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INTRODUCTION 


The research effort on the control of large space structures at the University 
of Houston has concentrated on the mathematical theory of finite-element models; 
identification of the mass, damping, and stiffness matrix; assignment of damping to 
structures; and decoupling of structure dynamics. The objective of the work has been 
and will continue to be the development of efficient numerical algorithms for analy- 
sis, control, and identification of large space structures. The major consideration 
in the development of the algorithms has been the large number of equations that must 
be handled by the algorithm as well as sensitivity of the algorithms to numerical 
errors. 

The finite-element model that has been used in the linear second-order matrix 
differential equation 


MdV tl + Cdx(t! + ^ (t) . £(t) (1) 

dt 2 dt 


where MeR mxm is the mass^matrix, CeR m><rn is the damping matrix, KeR mXm mx :j.s t * ie 
stiffness matrix, x(t)sR is the node displacement vector, and f(t)eR is the 

forcing function vector. 

The Laplace transform of equation (1) gives the matrix equation 


[Ms 2 +Cs+K]x(s) = B(s) 


( 2 ) 


where B(s) contains the initial condition information as well as the forcing func- 
tion. If s is replaced in equation (2) by X, equation (2) then takes the form of 
a lambda matrix: 


A(A)X(A) = B(A) 


( 3 ) 


If it is assumed that the initial conditions are zero and no forcing function is 
present, then 


A(A)X(A) = [MA 2 +CA+K]X(A) 


( 4 ) 
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is the homogeneous equation that will be of interest. Hie latent roots of 

A(X) are given by 


det A(A) = det[MA 2 +CA+K] = 0 


( 5 ) 


and the latent vectors y^ are obtained from 

A(X 1 )y i = [ MX 2 +C A ± +K] y ± = 0 


provided that all latent roots are distinct, of multiplicity one. 

Lancaster (ref. 1) and Dennis et al. (ref. 2), as well as others, have published 
material on lambda matrices. One of the tasks during the past research period has 
been to consider their work as well as extensions of the algebraic theory of lambda 
matrices for the control of structures. 


ALGEBRAIC THEORY OF LAMBDA MATRICES 

A comprehensive treatment of the algebraic theory of lambda matrices cannot be 
presented in this short paper. Only the essentials necessary to understand the damp- 
ing assignment problem will be given. 

Consider the lambda matrix 

A(A) = Q[IA 2 +CA+K]Q T = QA(A)Q T (6) 


T 

where Q is the Cholesky matrix of the decomposition of M = QQ . The normalization 
process in equation (6) will retain the symmetry of C and K. The lambda matrix 
A(\) will be considered in the following work and it will be assumed that the latent 
roots of A(\) appear in conjugate pairs as well as being distinct. This restric- 
tion is not necessary but is made to simplify the analysis. 

It can be shown that the matrix 


A = 
c 



(7) 
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has eigenvalues 
eigenvectors of A c 


that are equal to the 
contain the right and 


latent roots of A( A ) . 
left latent vectors of 


Hie right and left 
A(A) with 


and 



"ci 


(A i I+C)z i 



( 8 ) 


(9) 


where y . and z . are the right and left eigenvectors, respectively. An eigen- 

Cl Cl 

projector, P^, will be defined as the matrix 


P. 

l 


y .z . 
ci ci 

T 

Wei 


( 10 ) 


where the matrix A • defined by 


A . = P. A = A P. (11) 

Cl 1C Cl 

will have the same eigenvectors as A c but all eigenvalues will be zero except for 
which will be the same as in A c » The eigenprojectors are the matrix residues of 
the partial fraction expansion of [A c (\)] - ' 1 ' with 


[A c (X) ] 1 = 


2 m 

l 


i=l 



( 12 ) 


where 


- u-y'vvr 1 


(13) 
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A latent projector will 
tion expansion of [A(\)] _ ^ 


be defined as the matrix residues of the partial frac- 
where 


with 


[A(X)]" 1 


2m P ± 




(14) 


(15) 


The latent projectors are also given by 


A 




y i Z i 


dA(A ± ) 

dT~ 


(16) 


where y^ and are the latent vectors of A( X ) . 

To complete the limited discussion of the algebraic theory, the relationship 

that exists between the eigenprojectors and latent projectors must be given. By 

inverting A (\) it can be shown that 
c 


P 


i 



- 

A^+C I 
-K A ± I 


(17) 


with the past definitions holding. 


ASSIGNMENT OF DAMPING 

An algorithm to assign damping to a particular undamped mode has been developed, 
although the theory has not been fully explored at this time. Considerable work 
remains to be carried out to develop and test a comprehensive algorithm. Therefore 
the material in this section is only preliminary. 

Consider the undamped structure with the associated lambda matrix 

A(A) = IA 2 +K (18) 
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If it is assumed that K is symmetric and positive definite then all latent roots of 
A(\) lie along the jw ^xis. Since K is positive definite all of its eigenvalues 
will be real with \ = w £ , and the eigenprojectors of K are given by 


i m p -t 

(AI-K)" 1 = l - zr ^~ 
1=1 (A-A 1 ) 


( 19 ) 


where it can be shown that is given by 


— ~ 

P. = P. A.+P. 
i x i i 


( 20 ) 


which will be defined as the augmented latent projector. It is not difficult to show 
that the undamped A Q is 



( 21 ) 


but since the undamped matrix is 



( 22 ) 


then 


m 


l 

i=l 


A *9 

(p.a;+p,a/) = o 

l i i x 


( 23 ) 


m 


l 

i-1 


A 

< p iWi> 


= I 


( 24 ) 
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The damped case has a similar form with the exception that the upper diagonal block 
is modified with 



(25) 


Suppose now that the undamped structure is to have damping added to the 

it 

jth mode (i.e., Xj and \j) and all other modes are to remain undamped. Using 
equations (21) and (22), A c becomes 



/\ y\-k 


-<Yj +p iV k 


P.xf+P.x/ 

j 3 i j 


which can be rewritten as 


0 

m 

m 


• “ 

I P i 

i=l 1 

+ 

0 P. 

3 

"I 

i=l 1 

0 


-P.K -C. 
3 3 


- 


- 


(26) 


(27) 
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and also in the form 


cd 


= A 


cd 


2m 

I 

i=l 

±*3 


p i + 


A ,P. 
cd j 


+ 

A +A, 


(28) 


The matrix A + denotes the matrix constructed from the summation and Aj is the 
complement to A + . 

The last step in the algorithm is to recognize that A + is orthogonal to Aj ; 

thus 


r m _ - 


” — 


• • 

lf^' 

o 


0 p. 


0 0 

i=l 1 


J 








m 





-y p ,k o 


-P.K -C. 


0 0 

L 1 J 


_ J 3 


m m 


i=l 


(29) 


m 


Finally, C- is the matrix that is orthogonal to P,. The eigenprojector 

3 i=l 1 

Pj must be orthogonal to all as the set of eigenprojectors have the orthogonal 


property 


P_,P. = 0 
i J 


i ^ j 


(30) 


The Cj matrices are therefore nothing more than a.P . where a. is a scaling 
constant. The trace of Cj is 3 3 3 


m 

trace (C ) = £ c = 2a (31) 

J J 
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where trace (P . ) 
of A(\) with ^ 


1. The scaling constants are then used to place the latent roots 


'j 

2 


to, 


2 

(32) 

2 2 

t» . - CT. 

nj j 

(33) 


An example will now be given to illustrate the algorithm. Let C 
select K as 


0 and 


K = 


300 -200 

-200 350 


which has eigenvalues ±j 11.1105 and ±j22.9468. The augmented projectors are 


_ 0.562017 0.496139 

P = 

A 0.496139 0.437983 

and 

- 

_ 0.437983 -0.496139 

P 2 = 

-0.496139 0.562017 


The matrices 


C 


1 


1 

0.8827824 


0.8827824 

0.7793048 
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and 


C 


2 


1 -1.1327785 

-1.1327785 1.283187 


are orthogonal to and p^ . Note that the columns in and C 2 can be scaled 

by constants and remain orthogonal to P 1 and 

The eigenvalues of 


A = 
c 


-K ~C . 


are given in the computed results below. 


A(I,J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 

0 

0 

0 


0 

1 

0 

0 


EIGENVALUES 

1 

2 

3 

4 


REAL 

0 

0 

3.04403E-22 

3.04403E-22 


IMAGINARY 

22.9468 

-22.9468 

11.1105 

-11.1105 


No damping 


A{ I, J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 

0 

-1 

-.882783 


0 

1 

-.882783 

-.779305 


EIGENVALUES 


REAL 


IMAGINARY 


1 

2 

3 

4 


3. 99893E-08 
3.99893E-08 

- .889652 

- .889652 


22.9468 

-22.9468 

11.0748 

-11.0748 
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Damping in lowest mode 



A(I,J) MATRIX 


0 

0 

-300 

200 


0 10 
0 0 1 

200 -1 1.13278 

-350 1.13278 -1.28319 


EIGENVALUES 

1 

2 

3 

4 


REAL 

-1.14159 
-1.14159 
2 . 84795E-08 
2.84795E-08 


IMAGINARY 

22.9184 

-22.9184 

11.1105 

-11.1105 


Damping in highest mode 


A( I , J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 

0 

-1 

.125 


0 

1 

.125 

-1.03125 


EIGENVALUES REAL 


IMAGINARY 


1 

2 

3 

4 


-.570798 

-.570798 

-.444827 

-.444827 


22.9399 

-22.9399 

11.1019 

-11.1019 


Damping in both modes 


The undamped system matrix and its eigenvalues are given first. The lowest mode 
was then damped with the system matrix and its eigenvalue given. The third test was 
to include damping for the highest mode where C = C 2 . The last test was to combine 
the damping for the two modes with 


' 2 < C 1 +C 2 ) ‘ 


“l P l + “2 P 2 


0. 889655P 1 + 1.141595P 


The system matrix and its eigenvalues are then given. It should be noted that the 
new values of Oj of the two modes are 


0.889652 
G 1 = 2 


0.444827 
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and 


„ _ 1.141595 

°2 2 


.570797 


The algorithm allows damping to be assigned to all modes and to the prediction of the 
location of the system eigenvalues. 

The major problem that remains to be resolved is that of constructing the C 
matrix, which is symmetric and positive definite as well as realizable. It will be 
assumed that c and ic are tridiagonal with the off-diagonal elements having a 
smaller magnitude than the diagonal elements. Hie off-diagonal elements must be 
negative with positive diagonal terms. 


IDENTIFICATION OF M, C, AND K 

The quadrature algorithm for identifying the mass, damping, and stiffness 
matrices of a structure is still under study. Hie algorithm performs well with 
simulated data, but attempts to utilized data from the beam experiment at NASA 
Langley Research Center have not produced useable results. It is believed that the 
conditions on the beam during the data collection may not have satisfied the require- 
ments of the algorithm. 

The quadrature identification, as well as other algorithms, will receive major 
attention in the future. One of the major problems in developing an efficient 
algorithm is the availability of test data, either by computer simulation or from a 
test bed. Hie Langley Research Center experimental beam is a suitable test bed, but 
planning and running test data are costly and time consuming. Plans for the future 
include construction of an electronic analog test bed that will be low cost and will 
provide flexibility in the types of structures that can be simulated. Hie test bed 
wiilll be used primarily in the development of the identification algorithm, but it can 
also be used for developing control algorithms. 


SIMULATION OF STRUCTURES 

Some of the preliminary work on designing a structures test bed at the 
University of Houston has been completed. Hie type of simulator considered is an 
electronic analog that will be constructed from low-cost operational amplifiers, 
resistors, and capacitors. Data acquisition will be handled by a PDP 11/70 computer 
for signal processing. Hie PDP 11/70 is equipped with analog- to-digital and digital- 
to-analog converters and has adequate storage to perform on-line data collection and 
processing. 


CONCLUSIONS 

Some of the work that has been supported under grant from NASA Langley Research 
Center has been described in this report. Preliminary details on the assignment of 
damping are covered, as well as some information on identification and simulation. 
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The identification work has not been described since the algorithm i3 given in refer- 
ence 3. Hie simulation facility is still in the planning stage and a decision to 
build such a test bed has not been made. Further descriptions of this research are 
given in references 4 and 5. 
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